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The apparent accelerating expansion of the Universe, determined from observations of distant su- 
pernovae, and often taken to imply the existence of dark energy, may alternatively be explained by 
the effects of a giant underdense void if we relax the assumption of homogeneity on large scales. Re- 
cent studies have made use of the spherically- symmetric, radially-inhomogeneous Lemaitre-Tolman- 
Bondi (LTB) models to derive strong constraints on this scenario, particularly from observations 
of the kinematic Sunyaev-Zel'dovich (kSZ) effect which is sensitive to large scale inhomogeneity. 
However, most of these previous studies explicitly set the LTB 'bang time' function to be constant, 
neglecting an important freedom of the general solutions. Here we examine these models in full 
generality by relaxing this assumption. We find that although the extra freedom allowed by varying 
the bang time is sufficient to account for some observables individually, it is not enough to simulta- 
neously explain the supernovae observations, the small-angle CMB, the local Hubble rate, and the 
kSZ effect. This set of observables is strongly constraining, and effectively rules out simple LTB 
models as an explanation of dark energy. 



I. INTRODUCTION 

Observations of distant type-la supernovae are often 
taken to imply that the Universe has entered a phase of 
accelerating expansion, and may therefore contain 'dark 
energy' [TJ|2]. Such a conclusion, however, cannot be in- 
ferred from the supernova data alone - a model of the 
Universe is also required. At present, the simplest and 
most widely applied cosmological models are based on 
the Friedmann-Lemaitre-Robertson- Walker (FLRW) so- 
lutions of general relativity. These solutions are highly 
symmetric, and determining their validity as models of 
the real Universe is of critical importance for determin- 
ing the veracity of the claims involving dark energy. It is 
toward these ends that this study is aimed: Are spatially 
homogeneous and isotropic models with dark energy the 
only ones capable of accounting for the recent cosmolog- 
ical observations that appear to imply acceleration? 

The inference of acceleration is of profound conse- 
quence, not just for cosmology and gravitational physics, 
but also for particle and high energy physics. An acceler- 
ating Universe has an entirely different causal structure 
from one that is decelerating, with the vacuum itself tak- 
ing on a non-zero energy density and becoming thermal. 
Beyond this there is the 'cosmological constant problem', 
that contributions from the zero-point energy of quan- 
tum fields, and any bare cosmological term in Einstein's 
equations, must cancel up to 1 part in 10 120 [3 j. Such in- 
credible fine-tuning is widely believed to signify nothing 
less than a crisis in modern physics, and so the task of 
verifying the assumptions that go into our cosmological 
models becomes one of the utmost importance. 

Here we focus on the problem of radial inhomogeneity, 
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as modeled by the Lemaitre-Tolman-Bondi (LTB) solu- 
tions of general relativity [4]-[6]. These are the general 
spherically-symmetric solutions of Einstein's equations 
with dust. They are widely known to have more than 
enough freedom to account for the supernova observa- 
tions without recourse to dark energy [7], and are often 
referred to as 'void models' in the literature (but see also 
[8]). The relevant question is then whether or not these 
models are compatible with other observational probes 
of cosmology. 

This question has been addressed by various authors 
in a number of different contexts [9-21 . Most of these 
studies have, however, limited themselves to the special 
case of space-times that have a spatially homogeneous 
energy density at early times. This is achieved by con- 
sidering only those models that have a constant 'bang 
time'. In this case it is known that while the small angle 
CMB generated by a power-law spectrum of initial fluc- 
tuations can be easily reproduced within void models, an 
unacceptably low value of Ho is required to do so [20] [21] . 
However, it is also known that this problem can be alle- 
viated by allowing for general radial inhomogeneity, with 
non-constant bang time [20]. Here we address the prob- 
lem of whether or not other cosmological observables are 
also consistent with models that allow for this additional 
freedom, as well as further investigating the parameter 
space of solutions that fit the small- angle CMB. 

We will be interested in particular in the kinematic 
Sunyaev-Zel'dovich (kSZ) effect [22]. This effect is due 
to the inverse Compton scattering of photons from the 
Cosmic Microwave Background (CMB) off of electrons in 
distant clusters of galaxies. The rescattered light can be 
collected by observers who are looking at the cluster. If in 
the rest frame of the electrons the CMB has a non-zero 
dipole moment (in the direction of the observer) then 
the reflected light that the observer sees has its spectrum 
shifted. Such shifts are expected to be observable by 
upcoming experiments, and although they have yet to 
be directly detected [62], constraints have already been 
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placed on the allowed magnitude of this effect 

The kSZ effect is a particularly powerful probe of inho- 
mogeneity as it allows us to make observations not only 
along the null cone that is the boundary of our causal 
past, but also along null curves that go inside this cone. 
The distant galaxy clusters essentially act as mirrors, re- 
flecting light from the last scattering surface that would 
otherwise be unobservable to us. This extra information 
is above and beyond that which is available from the 
usual observations of distance measures, expansion rates 
and number counts, and so it is of great potential signif- 
icance as a cosmological probe. The power of the kSZ 
effect in this context appears to have first been pointed 
out by Goodman in [28] , although the first application of 
it to models that attempt to account for dark energy was 
performed by Garcia-Bellido and Haugb0lle [29]. These 
authors considered models with constant bang time only. 
We build on their work by allowing for a radially depen- 
dent bang time. 

To make progress it will be necessary to make a number 
of assumptions, which to avoid confusion we will state 
here. We assume the following: 

• That there is perfect spherical symmetry, with our- 
selves at the center of symmetry. 

• That the formation of the last scattering surface 
proceeds as in FLRW cosmology. 

• That there is a constant ratio of photons to baryons 
in the early universe. 

• That the spectrum of initial fluctuations is a power 
law in wave number, k. 

• That the energy density and all functions in the 
metric have smooth profiles. 

The first of these is inherent in the problem we have cho- 
sen to address. General perturbations to this exact set 
of symmetries have been considered in [30 and [53 , and 
the effects of being off-center have been considered in [19] . 
The second of these points is made for convenience. To 
date, we are unaware of any rigorous calculation involving 
the formation of the last scattering surface in inhomoge- 
neous space- times. The effect of allowing for an inhomo- 
geneous photon-to-baryon ratio has been considered in 
[3T] , and the related question of changing the position of 
the last scattering surface, while keeping the bang time 
constant, has been addressed in [18]. The effect of al- 
lowing a kink in the spectrum of initial fluctuations has 
been considered in [32]. We will not consider these free- 
doms further here, but note that a constant bang time 
is an assumption that would be added to similar lists in 
most other papers. For details of the effects of relaxing 
these assumptions, we refer the reader to the papers cited 
above. 

In Section [II] we present the LTB solutions, and dis- 
cuss how distance measures and redshifts are calculated 
within them. We then discuss the effects of the two radial 



degrees of freedom in these solutions, one of which is the 



bang time. In Section III we discuss some of the cosmo- 
logical probes that can be applied to these models, with 
particular reference to the kSZ effect. We also discuss 
why these observations are problematic for LTB mod- 
els with constant bang time. In Section [IV] we present 



our results, which include a detailed investigation of the 
effect of a radially varying bang time on CMB and Ho 
observations, as well as the kSZ effect. We show that, 
despite the additional freedom in the bang time, there is 
a combination of key observables that cannot be fitted 
simultaneously. We conclude in Section [V] that this ef- 
fectively rules out void models as an explanation of dark 
energy, unless one is prepared to discard one or more of 
the assumptions that we have listed above. 



II. THE MODEL 

In order to model general radial inhomogeneity we will 
use the Lemaitre-Tolman-Bondi (LTB) solutions of gen- 
eral relativity. These are given by the line-element [3H6] 



ds 2 = dt 2 



1 — k(r)r' 



; dr 2 -a 2 (t,r)r 2 dn 2 , (1) 



where ci2 = (air/, and where a\ must satisfy 



SttG m(r) k(r) 



(2) 



The functions k(r) and m(r) are arbitrary functions of 
the radial coordinate, and primes and over-dots denote 
partial derivatives with respect to r and t, respectively. 
These solutions are exact, and are the general spheri- 
cally symmetric dust-only solutions to Einstein's equa- 
tions. They admit a three dimensional group of Killing 
vectors that act transitively on the surfaces of constant 
r and t, and are spatially isotropic about the origin only. 

Solutions to Eq. ([2| are of the form a± — a\ (r, t — 
£b(0)> which introduces a third arbitrary function of r. 
This gives a total of three free functions: /c(r), m(r), 
and tB{r). We refer to these quantities as the spatial 
curvature, gravitational mass density (distinct from the 
local energy density) and bang time, respectively. In the 
limit of homogeneity they are all constant. It can also be 
seen that one can perform a coordinate transformation 
r — » f(r) that preserves the form of the metric in Eq. ([!]). 
This freedom can be used to set m =constant, without 
loss of generality (assuming (mr 3 )' is always positive). 
This leaves us with the general solution in terms of the 
spatial curvature, fc(r), and bang time function, ts(r), 
only. Analytic parametric solutions to Eq. ^ are known, 
and can be found in [33] . 

We can now use these solutions as cosmological models 
that exhibit an arbitrary amount of radial inhomogeneity 
by supposing ourselves to be observers at the center of 
symmetry. Such models are known to be able to produce 
excellent fits to the supernova data without requiring any 
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dark energy, and often result in the observer being at the 
center of gigaparsec-scale underdensity, or 'void'. This 
is possible due to both temporal and spatial variations 
in the geometry of the space-time that are experienced 
by photons as they travel through the void. Such calcu- 
lations require knowledge of redshifts and distance mea- 
sures in this space-time, which we will now consider. 

Let us first define two different Hubble rates: a trans- 
verse one, Hi = di/ai, and a radial one, H2 = 0,2/0,2. In 
the limit of homogeneity these two quantities are identi- 
cal, but differ, in general, in inhomogeneous space-times. 
The redshift of photons traveling along radial geodesies 
can then be calculated by integrating the radial Hubble 
rate as follows: 



1 + z = exp 



H 2 (t,r(t))dt 



}■ 



(3) 



where r = r(t) is a solution of the radial geodesic equa- 
tion, and t e and t Q are the time at which the photon was 
emitted and observed, respectively. Note that the rela- 
tion (1 + z) oc 1/ai no longer holds, in general. We can 
also calculate the angular diameter distance to objects at 
redshift z using 



d A (z) = r(z) a 1 (r(z),t(z)), 



(4) 



where t(z) is calculated by inverting z = z(t) from Eq. 
([3]), and r(z) — r(t(z)) is found using the radial null 
geodesic equation. Luminosity distances are then given 
by Etherington's reciprocity theorem [34] 



d L (z) = (l + z) 2 d A (z), 



(5) 



which is true in any space-time. The local energy density 
is given by p = (mr 3 )' / '3a20 2 r 2 . 

In this paper we will often choose to parametrize the 
two functions k(r) and ts(r) as Gaussian curves, with 

k(r) = A fc exp(-r 2 /A^) + /c oc (6) 
t B (r) = A tB eM-r 2 /\ 2 tB )f(r), (7) 

where Ak, A tB , A&, X tB and k^ are constants, and the 
factor f(r) = exp(— r 10 /A^) is included to attenuate 
the bang time profile at large r. This is done so that 
wide profiles can be used, while limiting the effect of 
early inhomogeneity on the central observer's last scat- 
tering surface, and is discussed further in Section |IIB| 
The profiles above are defined by their amplitudes, A^ 
and widths, A^. A further parameter k^ defines the 
asymptotic spatial curvature, outside the void. The 
timescale of the model is set by choosing a local Hub- 
ble parameter Hq = Hi\ r= o = #2^=0 at time t G , where 
(to — £#)|r=o is the age of the Universe along the world- 
line of an observer at r = 0. A rescaling can be used to 
set ai(0,t o ) = a 2 (0,t o ) = 1. 

The amplitudes in Eq. (|6| and ^ can be expressed in 
a more familiar form as fractions of the total density at 



the origin today. For this purpose let us define 

-A k /H 2 
~k QO / Hq 



&k 2 



^ m. — 



8nG , , 
3M m(r) ' 



(8) 
(9) 
(10) 

(11) 



such that fifc + Q m = 1. Furthermore, restricting our- 
selves to Vt m > means that we consider only < 1. 

Let us now briefly consider the consequences of fluctua- 
tions in k(r) and ^(r). The former of these is the analog 
of the spatial curvature in FLRW solutions, which dom- 
inates the dynamical evolution of the Universe at late 
times. The latter changes the location of the initial sin- 
gularity, and so can be thought of modifying the early 
stages of the Universe's history. This interpretation is 
supported by treating the LTB geometry as a fluctuation 
about an FLRW solution. In this case the fluctuations in 
k(r) can be mapped into growing modes, while fluctua- 
tions in ts(r) are mapped onto decaying modes [35] . 



A. An Inhomogeneous Late Universe 

Let us first consider k(r). It can be seen fromEq. Q 
that k < gives a positive contribution to the expansion 
of ai, while k > gives a negative contribution. This 
is the behavior we are familiar with from the Friedmann 
equation of FLRW cosmology. Unlike the homogeneous 
FLRW solutions, however, the expansion of 02 does not 
always get a positive contribution from k < 0. On the 
contrary, in regions of k < the expansion of 02 can 
slow, and recollapse can occur. This behavior is well 
known, and can lead to the formation of a 'shell crossing 
singularity' when the collapsing region reaches 02 = 0. 

One can avoid shell crossing singularities in a region 
by satisfying the Hellaby-Lake conditions [36] . These de- 
pend on the sign of fc(r), and for k < may be written 
using our notation as [37 



(mr 3 Y > 0, t' B < 0, and (kr 2 )' < 0, 



(12) 



while for k > the last of these should be replaced by 



log 



3* ; g|fc| 
87rGra 



> 0. 



(13) 



These conditions guarantee that 02 > 0, so that shell 
crossing singularities cannot occur [63]. They are, how- 
ever, very restrictive, and most applications of the LTB 
solutions to cosmology simply avoid the issue by mak- 
ing sure that shell crossings only happen in the distant 
future. They can then be considered as a breakdown 
of the model at some future time, after which a more so- 
phisticated solution including pressure would be required 
to avoid the formation of singularities. The existence of 
pressure is expected to prevent the complete collapse of 
matter, and a large over density of collapsed structures is 
thus expected to form instead. 
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B. An Inhomogeneous Big Bang 

Let us now consider the consequences of fluctuations 
in tsir). If tsir) is not constant, then the 'age of the 
universe' differs from place to place. This is a significant 
departure from the standard picture of the big bang, and 
may initially seem odd. Certainly, there have been a 
number of objections to allowing inhomogeneous bang 
times in the literature, with the result that to date most 
studies of void models have expressly set tsir) =constant 
a fortiori. In this section we will discuss the physical sig- 
nificance of an inhomogeneous big bang, and argue that 
it is reasonable to consider models with such a feature. 

As mentioned above, fluctuations in £^(0 correspond 
to decaying modes when the space-time is approximated 
as a perturbed FLRW solution. A non-constant bang 
time therefore corresponds to an inhomogeneous early 
universe, and as one goes further back in time the size 
of the consequent inhomogeneity generally increases. As 
with the case of fluctuations in &(r), there exist points be- 
yond which the scale factor is contracting rather than 
expanding, and in cases where the Hellaby-Lake condi- 
tions are violated, shell crossing singularities can occur. 
In the case of fluctuations in ts(r), however, this behav- 
ior occurs at very early times rather than at very late 
times. 

In Fig. [T] we illustrate the existence of surfaces of 
ci2 = in cases with t' B > 0, and regions with H2 < 
in cases with t' B < 0. The former of these are the early 
universe analogue of the shell crossing singularities we 
described in the previous section, and in this case we con- 
sider the singular surface with ci2 = to be our initial 
hypersurface. The latter case corresponds to regions that 
have started to collapse, but reach the singular surface 
t = tB before shell crossings occur. Contraction of this 
type should be expected to cause blueshifts when looking 
along exactly radial geodesies [38 , and again this behav- 
ior has an analog in the inhomogeneities that form at late 
times [9j [17] . Such blueshifts would have profound effects 
if they were allowed to occur between the last scattering 
surface and an observer (for example, a distance-redshift 
relation r(z) that is not monotonic could occur). 

Celerier et al. [8 found that small variations in the 
bang time, of order hundreds of years, lead to tempera- 
ture anisotropics in the CMB of (9(10 -6 ), which are cur- 
rently only marginally too small to be observed. Larger 
bang time variations would have a stronger observational 
signal, but this need not be an issue if the region of vary- 
ing bang time occurs far inside the void, away from the 
surface of last scattering that we see directly. Neverthe- 
less, observers elsewhere in the space- time would see con- 
siderable anisotropics in their CMB sky, and this could 
be observable in the kSZ effect we see from CMB photons 
that rescatter off their cluster. The kSZ effect therefore 
has the potential to provide powerful constraints on the 
inhomogeneities caused by bang time fluctuations. 

Of course, LTB models are dust-only solutions of Ein- 
stein's equations, and lose their validity when radiation 
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FIG. 1: Upper panel: If t' B > then a surface exists with 
a 2 — 0, corresponding to the occurrence of a shell crossing 
singularity (dotted line). Beyond this, the formal solution 
for the energy density gives negative values (hatched area). 
Lower panel: If t' B < 0, regions with H2 < form (hatched 
area). The solid lines correspond to the surfaces t — ts- 

becomes important at early times [8] [51] . The intro- 
duction of a radiation fluid into inhomogeneous solutions 
complicates matters considerably, and we do not attempt 
to include the gravitational effects of radiation in our 
models here. However, in the same way that one would 
expect pressure to prevent the formation of shell crossing 
singularities at late times, one could reasonably speculate 
on a similar mechanism occurring at early times. In the 
present study we concentrate on the matter dominated 
phase of the Universe's history, which should be suffi- 
cient to model the Universe from last scattering to the 
present time. This is well modeled by the LTB solutions. 
We leave the consideration of the gravitational effects of 
inhomogeneous radiation fields to other studies. For fur- 
ther details of this the reader is referred to [5TJ [60] . 

An obvious concern with models of this type is that 
they are difficult to reconcile with early universe infla- 
tion. This is true with models that have ts =constant, 
as well as models with an inhomogeneous big bang. One 
must then either discard inflation for the time being, or 
attempt to construct inflationary scenarios that result in 
occasional large inhomogeneities (see, e.g., [39. and [40]). 
Here we address the problem of what can be said about 
the geometry of the Universe directly from observations, 
rather than imposing requirements from theories of the 
very early universe. 



III. OBSERVATIONAL PROBES, AND 
RESULTS WITH CONSTANT BANG TIME 

Void models that reproduce the observed supernova 
distance modulus curve have proved relatively easy to 
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construct, and little more than a moderately deep un- 
derdensity with a comoving width of order a gigaparsec 
is required to obtain a satisfactory fit. Indeed, the ease 
with which the supernova observations can be reproduced 
has been one of the principle factors motivating interest 
in these models. 

The introduction of such a large inhomogeneity, how- 
ever, can hardly be expected to leave predictions for other 
cosmological observables unchanged, and so there have 
been a number of attempts to make detailed tests of 
voids using multiple data sets [9] [T3l l2T] . While thor- 
ough, these previous studies have limited themselves to 
the case of voids with constant bang times. In this sec- 
tion we summarize the constraints that can be imposed 
on such void models from observations of supernovae, the 
CMB, the local Hubble parameter, and the kSZ effect. In 
particular, we draw attention to the difficulty that voids 
with constant bang time have in fitting the CMB and Hq 
simultaneously, and discuss the power of the kSZ effect 
as a test of large-scale inhomogeneity in these models 

nana Esi . 



In Section IV we will proceed to consider more general 
voids with varying bang times. This additional freedom 
allows some of the constraints on the specific observables 
discussed in this section to be weakened significantly (al- 
though a combined fit to all data sets remains elusive). 



A. Supernovae 

As noted above, fitting the supernova data is a rela- 
tively simple matter, and void models can be constructed 
that fit any given (Il(z) curve exactly [41]. One should be 
aware, however, that reproducing the precise effects of A 
at low z requires an energy density profile that is 'cusped' 
at the center p~2j [42] . Generic smooth profiles produce 
qualitatively different behavior, due to the Milne-like ge- 
ometry near the origin, but can still be shown to be con- 
sistent with current data sets [T2| 142]. 



In general, voids and FLRW models with the same lo- 
cal geometry have different angular diameter distances 
to the last scattering surface, resulting in a relative shift 
in their observed CMB power spectra. Now, the distance 
to last scattering can be adjusted by changing the width 
and depth of the void, but this typically produces rela- 
tively small shifts that are not enough to bring the peak 
positions of the CMB power spectrum in line with cur- 
rent observational constraints [20]. Changing the curva- 
ture of the FLRW region near the last scattering surface, 
however, produces much larger effects [20], and good fits 
to the small-scale CMB power spectrum can be found 
for void models that have positive asymptotic curvature. 
Such models, however, require an anomalously low local 
Hubble rate (Hq < 50 km s _1 Mpc -1 ) in order to keep 
the expansion rate at last scattering low enough to be 
consistent with the data [20j [21]. This is strongly in- 
consistent with the observed value of Ho = 73.8 ± 2.4 
km s -1 Mpc -1 recently found in [43 . The CMB+# by 
themselves are therefore sufficient to effectively rule out 
simple void models with constant bang time. 

One can attempt to avoid this conclusion by violating 
the assumptions that we set out in Section [I] In particu- 
lar, models with inhomogeneous last scattering surfaces 
have been considered in [14] and [31], and a non-power 
law spectrum of initial fluctuations has been considered 
in [32] . If one is prepared to consider such additional free- 
doms then the CMB+i^o constraints can be considerably 
weakened. 

In Section |IV| we will consider the consequences of 
CMB observations in general void models, where the 
bang time is allowed to vary. It has already been shown 
in [20] that the available constraints from the CMB-hHo 
can be considerably weakened in this case, without any 
need to violate the assumptions introduced in Section [I] 
In Section [IV] we will quantify this result, finding best fit 
models and confidence regions in parameter space. We 
will show that the best fit models have bang time fluctu- 
ations of order a billion years. 



B. The CMB and H 

If the last scattering surface we observe is located in a 
region of the Universe that is homogeneous and isotropic 
enough to be modeled as being approximately FLRW 
then we can use standard techniques to calculate the 
power spectrum of fluctuations on that surface. The 
CMB that we measure on our sky then depends on the 
initial spectrum, which can be calculated using an ef- 
fective FLRW model, and the projection of fluctuations 
from the last scattering surface onto our sky. This projec- 
tion depends on the space-time geometry between us and 
the last scattering surface, and can be calculated from the 
angular diameter distance in Eq. Q. It is in this way 
that the CMB provides constraints on the geometry of 
the late Universe. 



C. The kSZ Effect 

A promising observable for testing large-scale homo- 
geneity is the kinematic Sunyaev-Zel'dovich (kSZ) effect. 
This effect occurs because galaxy clusters contain hot gas 
that can Compton scatter CMB photons, leading to a 
frequency-dependent temperature increment /decrement 
of the CMB along the line of sight to the cluster. Such 
scattering events cause two separate effects that can be 
distinguished in the reflected light. The first is due to 
transfer of thermal energy from the cluster gas to the pho- 
tons, and is known as the thermal Sunyaev-Zel'dovich ef- 
fect. The second is due to the dipole, AT/T, of the CMB 
radiation on the sky of an observer comoving with the re- 
flecting cluster, and is known as the kinematic Sunyaev- 
Zel'dovich effect [22]. It produces a temperature change 
of AT/T in the reflected light, and is similar to the rel- 
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ativistic Doppler shift that one would experience by re- 
flecting a beam of light off of a moving mirror. 

As noted by Goodman [28], and explicitly calculated 
by Garcia-Bellido and Haugb0lle [29], observers who are 
off-center in a radially-inhomogeneous universe should 
expect to see a large AT/T in their CMB if they are 
comoving with the dust. This is because the distance- 
redshift relation becomes a function of direction in an in- 
homogeneous space- time, so that the surface of last scat- 
tering will appear to be at different distances/redshifts in 
different directions on the sky of an off-center observer. 
The result of a large kSZ effect then follows because most 
observers in the space-time see a highly anisotropic CMB. 
Of course, this is not the case in an FLRW universe, 
where one should anticipate a low kSZ signal due only 
to the peculiar motion of clusters. It is for this reason 
that the kSZ effect is expected to be a powerful probe of 
large-scale inhomogeneity. 

In a void model the dipole in the CMB is aligned in the 
radial direction due to spherical symmetry, and can be 
calculated from the relative velocity to an observer at the 
same point who would see an isotropic CMB. It is this 
dipole that can then, in principle, be measured using the 
kSZ effect. Now, the magnitude of dipole, AT/T, can be 
calculated for a given void model by finding the redshifts 
to last scattering when looking radially into and out of 
the void. The observer then sees an average temperature 
of T = \{T in + T ont ), where T in and T out are the temper- 
atures of CMB photons seen when looking into and out 
of the center of symmetry, respectively, and the relative 
velocity with respect to the CMB rest frame that causes 
this dipole is given by 



AT 



Zout 



%out 



(14) 



where Zi n and z out are the redshifts to last scattering in 
the directions toward and away from the center of the 
void, respectively. They can be calculated from Eq. ([3]), 
which is valid for off-center observers [44]. The kSZ effect 
can also be measured as a power spectrum using [15] fl6] 



AT(n) 



T 



kSZ 



z * dr 

/3(z)5 e (n,z)—dz, (15) 



where AT(h)/T is the CMB anisotropy seen by an off- 
center observer in a direction n, r is the optical depth 
along the line of sight, S e is the density contrast of elec- 
trons, and z* is the redshift to the last scattering surface. 

One should note that the calculation described above 
over-estimates AT/T in the reflected light because the 
anisotropy seen by off-center observers will not be purely 
dipolar, especially far from the center of the void [19]. 
The dipole contribution, however, is the dominant one at 
low z, and so we expect the prescription outlined above 
to be accurate enough for our current purposes. Obser- 
vations of individual clusters have yielded upper limits of 
AT/T < 2000 kms -1 [29 , and more recent observations 
from ACT and SPT have produced upper limits on the 



kSZ power spectrum at I = 3000 of 8 fiK 2 and 13 /iK 2 , re- 
spectively [26] [27] . This is consistent with the typical pe- 
culiar velocities expected in ACDM of ~400 kms -1 , but 
is strongly inconsistent with any large void with constant 
bang time that obeys the assumptions made in Section [T] 

An example f3(z) profile that an observer at the cen- 
ter of a large void with constant bang time could infer 
from observations of the kSZ effect is shown in Figure [2] 
Although this is only one example, the enormous magni- 
tude of the effect is a generic result for observers located 
at the center of such voids. This directly demonstrates 
the utility of the kSZ effect as a probe of inhomogeneity 
on large scales, and explains why current observational 
constraints on the kSZ effect by themselves are enough to 
rule out simple voids with constant bang time. In Section 
IV we consider the consequences of a varying bang time 
on observations of the kSZ effect, and show that there 
exist general giant void models that are consistent with 
constraints from current observations. 




FIG. 2: An example of the relative velocity with respect to 
the CMB rest frame, /3(z), that would be inferred by a central 
observer in an example giant void with constant bang time. 



D. Other Observables 

So far we have discussed the specific observables of 
supernovae, the CMB+i^o, and the kSZ effect, and de- 
scribed how the latter two of these provide constraints on 
simple void models with a constant bang time that are 
sufficient to effectively rule them out. We have chosen to 
discuss these particular observables as we consider them 
to be reasonably well defined in void models, easily calcu- 
lable, and very constraining. There are, of course, other 
observables that one could also consider. These include 
baryon acoustic oscillations (BAOs), galaxy correlation 
functions, and the integrated Sachs- Wolfe (ISW) effect, 
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to name just a few. 

To understand these observables in the FLRW cosmo- 
logical models one only needs linear perturbation the- 
ory about an FLRW background, which is very well un- 
derstood. The observables in question can therefore be 
straightforwardly calculated. Linear perturbation theory 
in LTB cosmology, however, is significantly more com- 
plicated. In [30 a gauge invariant formalism for gen- 
eral perturbations in spherically symmetric space-times 
is applied to these models, and it is shown that scalar, 
vector and tensor modes no longer decouple. This means 
that complicated effects can occur that are not present in 
FLRW cosmology. The full consequences of this behav- 
ior have yet to be understood, and so here we avoid the 
use of observables that rely on linear perturbation theory. 
This includes BAOs, galaxy correlation functions and the 
ISW effect. For treatments of some of these observables 
when t b =constant the reader is referred to [T7] . 

IV. RESULTS WITH VARYING BANG TIME 

In this section we examine the constraints that can be 
imposed on general void models in which the bang time 
is allowed to vary. This generalizes the previous results 
that were summarized in Section Hill 

As before, the observables we will use to constrain 
these models are the supernova distance moduli as func- 
tions of redshift, the CMB power spectrum on small 
scales, the local Hubble rate, and the kinematic Sunyaev- 
Zel'dovich effect. The specific data used for each ob- 
servable will be explained in the subsections that fol- 
low. We use the parametrized LTB models of Section 
|lH and a Metropolis-Hastings Markov chain Monte Carlo 
(MCMC) method to explore parameter space [64 . The 
likelihood function for each set of parameters is modeled 
as a chi-squared distribution, such that — 21og£ « x 2 , 
and 'goodness of fit' is quantified by comparing to a 
ACDM model with Vt K = 0.734 and h = 0.710 [45]. 

We first proceed by considering each observable in- 
dividually, and then go on to consider the constraints 
available from combinations of different observables. It 
is found that the additional freedom allowed by varying 
the bang time significantly weakens the constraints that 
each observable imposes by itself, but that the combined 
power of all observables is still enough to effectively rule 
out these models as a possible explanation of dark en- 
ergy. In particular, we show that neither the CMB+i^o 
observations, nor the upper bounds on the kSZ effect 
for individual clusters, have the ability to rule out these 
models by themselves, as is the case when the bang time 
is assumed to be constant. 



A. Super novae 

In the fits that follow we use the Union2 compilation 
of 557 supernovae, which extends out to z ~ 1.4 [46] . 



Other supernova data sets also exist, and void fitting 
procedures are known to exhibit some sensitivity to the 
data set that is chosen [12, 32 . We choose the Union2 
data as it is the most extensive catalog, and the most 
widely used in the literature. The absolute magnitude 
of the supernovae in the Union2 data set is an unknown 
parameter, and is therefore fitted to each model individu- 
ally as a nuisance parameter. We use the published errors 
in this data set, which includes an 'intrinsic error' that 
is added to minimize the reduced y 2 of ACDM. The full 
Union2 "covariance matrix with systematics" is used in 
performing all of the likelihood estimates that follow. 

As with ts =constant, there is no problem fitting the 
supernova data without dark energy. 

B. The CMB and H 

In FLRW cosmology, and for our current purposes, 
the CMB power spectrum can be efficiently specified on 
small scales with only three pieces of information [65] : (i) 
the acoustic horizon scale at decoupling, (ii) the acoustic 
scale at matter-radiation equality, and (iii) the projected 
scale of the CMB onto our sky. This information can be 
combined into three parameters in a number of different 
ways [3TJ [47H49] , but here we choose to specify it as the 
'shift parameter', 5, the Hubble rate at last scattering, 
H*, and the redshift of the last scattering surface, z*. 
The shift parameter is defined as S = g?a(2*)/^U(2*)> 
where is the angular diameter distance to the last 

scattering surface in a fiducial spatially flat FLRW model 
with Qm 1. This quantity corresponds to the change 
in scale of fluctuations on the sky that two observers in 
different space-times would see when looking at two iden- 
tical last scattering surfaces. 
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FIG. 3: Marginalized likelihood for the parameters S and 
if*, found using WMAP 7- year data and a modified version 
of CosmoMC [51]. Shaded regions show the 68%, 95%, and 
99.7% confidence regions. 
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FIG. 4: Marginalized likelihoods for a Gaussian void model with varying bang time. Shaded regions show the 68%, 95%, and 
99.7% confidence regions. Here we use Ho = lOO/i kms -1 Mpc -1 . See the text for other definitions. 



Here, for simplicity, we take z* = 1090, which is the 
redshift to the last scattering surface in ACDM and var- 
ious other models [48] [66]. It now remains to impose 
constraints on S and H*. To do this, we enforce the 
condition that the region of space in which the last scat- 
tering surface forms is well- approximated as being ho- 
mogeneous and isotropic, so that standard results from 
FLRW cosmology can be applied in dealing with all of 
the physics up until the formation of the last scattering 
surface. We can then use CosmoMC [50] to calculate H* 
and (Ia(z*) for an observer in a spatially flat and dust 
dominated FLRW universe looking at this surface. Us- 
ing the space-time geometry of our void models we can 
then calculate S and Hq for an observer at the center of 
the void looking at an identical last scattering surface, 
with identical Hubble rate at last scattering, and at an 
identical redshift. This is the procedure followed in [20] . 

We use the WMAP 7- year data [51], with a modified 
version of CosmoMC, to constrain our models. In this 
analysis we choose to only use data at £ > 100, as the 
\ow-£ power spectrum is dominated by the ISW effect 



(see Section HID for a brief discussion of this). This 



choice weakens the constraints that can be achieved on 
the scalar spectral index of the initial power spectrum, 
n s . Conservatively, we fix n s = 0.96 here [67]. The 
constraints that can then be imposed on S and H* are 
shown in Fig. [3] The best fit values are found to be S — 
0.875 and ~ 1.27 x 10 6 kms^Mpc -1 . These values 
are consistent with those found in [20 , 21 . Note that the 
CosmoMC CMB fits were not performed jointly with the 
void model MCMC; instead, they were run beforehand 
to get likelihoods for S and which we then used as 
priors for the void model MCMC. 

As discussed in Section 



IIIB[ it is possible to con- 



struct simple void models with =constant that sat- 
isfy the constraints on S displayed in Fig. [3] This can 
be achieved by simply changing the spatial curvature 
of the model at large z [20 j. The constraints on H*, 
however, are more difficult to satisfy. For simple Gaus- 
sian voids with ts =constant, under the assumptions 
described above, the WMAP 7- year data [51 and the 
Union2 supernova data set [46] are enough to show that 
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Hq < 40 kms _1 Mpc _1 is required, which is in strong dis- 
agreement with the value of Ho = 73.8±2.4 kms _1 Mpc _1 
found by Riess et al. [43 j. Allowing z* to vary can in- 
crease the upper bound on Ho by around 5 kms _1 Mpc _1 , 
and changing the precise functional form of k(r) can also 
marginally change Hq (see [9j[T7]). These are relatively 
small effects, however, and unless one is prepared to re- 
ject one or more of the assumptions given in Section [I] 
models with ts =constant remain strongly inconsistent 
with recent measurements of Hq. 

Allowing the bang time function to vary significantly 
improves the ability of void models to fit the CMB+iiZo 
data [20]. In Fig. [4] we show the likelihood plots 
for the parameters fi^ , A tB , X tB ,Hq), when con- 
strained with the WMAP 7- year data [51], the Union2 
data set [46], and the measurement of Ho — 73.8 ± 2.4 
kms _1 Mpc _1 [43 . Good fits to the data are obtained for 
models with a bang time fluctuation of width 8000 Mpc 
that makes the Universe about 800 million years older 
in the center than it is at large r. The curvature profile 
is narrower than this, with a width of 2500 Mpc, and a 
depth of Qk — 0-83 at the center. The preferred spatial 
curvature at large radii is only Qk,2 ~ +0.002. It can 
be seen that in this case the model is able to produce 
an acceptably large value of Hq, with a best-fit value 
of 73.6 kms _1 Mpc _1 . When compared to the best-fit 
ACDM model we find that ACDM is slightly preferred, 
with A% 2 =+4.5 for 560 degrees of freedom. Most of this 
difference is due to the void model having a poorer fit to 
the supernova data, even though it over-fits Ho and the 
CMB data. Voids with slightly more complicated spa- 
tial curvature profiles produce fits to the data that are 
at least as good as ACDM. 

In the best fitting models the curvature profile, fc(r), is 
largely responsible for shaping the void at low redshift, 
with z < 1. In this region the bang time gradient is 
small, and so has little effect. In the region 1 < z < 2 
the curvature profile then flattens out and the bang time 
gradient begins to change rapidly. This produces large 
fluctuations in H2(z) along our past null cone, such that 
H 2 can take lower values at large z. The low value of H* 
required at last scattering can then be simultaneously 
accommodated with a large value of Ho locally. The 
difference in profile widths therefore helps to explain how 
it is that a good fit to the data can be achieved. 



C. The kSZ Effect 

Let us now consider the kSZ effect in void models with 
varying bang times. The void-induced dipole, AT/T, can 
be calculated in an LTB model by following the procedure 
below: 

1. On every point on our past null cone, solve the 
radial null geodesic equation for light rays traveling 
both into and out of the center of the void. 

2. Calculate the redshift to the last scattering surface 



along these geodesies using Eq. (§. 



3. Calculate the dipole, AT/T, using Eq. This 
can be converted into an effective velocity using 
ft = AT/T. 

This procedure relies on knowing the location of the last 
scattering surface at different values of r. For models 
with a constant bang time, this surface occurs at a con- 
stant time, t = tLS- I n models with varying bang time, 
however, it will not occur as a hypersurface of constant 
t, as the presence of a bang time gradient changes the 
time evolution of the radial Hubble rate and density at 
a given r. We therefore approximate the location of the 
last scattering surface as a hypersurface of constant den- 
sity, p, rather than time, t. The precise location of the 
last scattering surface will turn out to be important, and 
we will discuss the consequences of altering its position 
as we proceed. 

We use the upper limits on ft = AT/T that have been 
measured from nine individual clusters [23-25 , as col- 
lected in [29 . These clusters span a redshift range of 
0.18 < z < 0.55. The data have asymmetric statistical 
errors, and are subject to large systematic errors of up to 
~ 750 kms -1 . For further explanation of the uncertain- 
ties in this data set the reader is referred to [29] . 

Void models with constant have already been shown 
to be inconsistent with even this limited data set, as we 



discussed in Section IIIC 29 . Within this class of mod- 



els, and subject to the assumptions outlined in Section 
[T] the best fitting voids are those that are either very 
shallow (flfa <^ 1) or ver y narrow (A& <C 1 Gpc), with 
the latter of these possibilities only working because it 
restricts the inhomogeneity to redshifts at which there 
are currently no data. We find that the extra freedom 
afforded by allowing the bang time to vary relaxes these 
tight constraints, and admits the possibility of allowing 
ft(z) to be small even in regions of the Universe that are 
strongly inhomogeneous, with as large as 0.85. 

Fig. [5] shows an example of a large void with varying 
bang time that produces small enough ft(z) to be com- 
patible with the data discussed above. A model with the 
same curvature profile, fc(r), but a constant bang time is 
also displayed. It can be seen that the additional free- 
dom allowed by the varying bang time has a considerable 
impact on ft(z). The energy density profile for this model 
is also displayed in the figure. One should note, however, 
that the functional form of the bang time fluctuation re- 
quired to produce this result is more complicated than 
the simple profile of Eq. [7[ Instead, a sum of two (mod- 
ified) Gaussian curves was used, of the form 

t B (r) = A, exp (- [(r - n) 2 /A? + (r - n) 8 /Af]) (16) 
+ A 2 exp (- [(r - r 2 ) 2 /A| + (r - r 2 f/\t]) , 

where Ai$, ^1,2, and Ai 5 2 are the amplitude, offset from 
the origin and width of the Gaussians respectively. De- 
spite the greatly increased freedom in this bang time 
profile, we were unable to find a model consistent with 
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FIG. 5: Upper panel: Velocity with respect to the frame 
where the CMB is isotropic, /3(z), for models with constant 
bang time (dashed blue line) and non-constant bang time 
(solid red line). Both models have the same spatial curva- 
ture, k(r). The data points are the upper limits for the nine 
clusters used in [29]. The solid red curve has a bang time 
fluctuation that is the sum of two Gaussians. Lower panel: 
Normalized density as a function of r on a hypersurface of 
constant t for the same two models. 



the kSZ and supernova data simultaneously; the best- 
fit model had AxIn = 4 4 - 1 witn ACDM. It is plausible 
that the situation could be improved by considering yet 
more complicated functional forms for k(r) and £b(t*), a 
possibility that we investigate in Section |IVE| For sim- 
plicity, we have assumed in Fig. [5] that the offset caused 
by systematic errors in the data is zero. In reality, the 
data points would likely move slightly toward the curve 
to which they are being fitted, in order to improve the 
likelihood. 

Finding models that agree with the upper limits on 
the statistical kSZ from ACT and SPT is a more difficult 
task. The kSZ power spectrum given by AT 2 , in Eq. 



(15), depends on an integral of f)(z) over redshift. Devi- 
ations from P = at any redshift therefore accumulate, 
potentially producing a large kSZ signal. It is possible 
that /3(z) could be made to change sign so that nega- 
tive contributions cancel the positive ones, but this would 
require a delicate balancing of the competing effects to 
satisfy the statistical and single-cluster kSZ data simulta- 
neously. One should bear in mind, however, that at large 
enough z, the effect of the void on the observed kSZ ef- 
fect will decrease, as the angle subtended by the void 
on the distant observer's sky decreases, and the power 
in the dipole term of the anisotropy is shifted to higher 
multipoles [T9] . 



In summary, we find that large void models with vary- 
ing bang times may have enough extra freedom available 
to alleviate the constraints that can currently be imposed 
from observations of the kSZ effect. As we discuss in the 
following section, however, it is unlikely that after doing 
this there will be enough remaining freedom to accom- 
modate any other observables. 



D. Combined Constraints 
(SN+CMB+# +kSZ) 

Let us now consider combining all of the observables 
we have discussed so far. These are the Union2 supernova 
data, the WMAP 7- year data, local measurements of Hq, 
and the kSZ effect. 

In Fig. [6] we show the observed value of AT/T as 
a function of z that a central observer would measure 
from the kSZ effect in the models found in Section [IV Bl 
These models have been shown to provide good fits to 
the supernova data, and the CMB and Ho data sets si- 
multaneously. In Fig. [7] we show this information as a 
function of H for redshifts z = 0.05, 0.10, 0.15, and 0.20. 
At all redshifts considered the distribution was bimodal, 
with some models having j3 ~ 1. Such incredibly high 
velocities are completely inconsistent with the data, and 
so here we show only the models with lower (3. It can be 
seen that the value of the kSZ signal that one would ob- 
serve from the center of these models is extremely large, 
even at low redshift. Such enormous kSZ signals are not 
compatible with the data displayed in Fig. [5j even with 
very large additional systematic uncertainties included. 

The principal reason for this large effect appears to be 
the large width of bang time fluctuation that is favored 
by the combination of supernova, CMB and Ho data sets 
(see Fig. [Z]). Because of this, observers at z > 0.1 look 
through regions in which the bang time gradient is large 
when they look through the void. As discussed in Section 
II B| these regions host shell crossings when t' B > 0. This 
pushes the surface of last scattering to much later times, 
and causes significant modifications to the redshift that 
this surface is seen at when looking through the center of 
the void. The redshift of the last scattering surface when 
looking away from the void experiences no such effect, 
as it is effectively fixed in position by the CMB data we 
see from the center. As a result, the values of Zj„ and 



z out in Eq. 14 differ significantly, and the value of AT/T 
is therefore even larger than in the constant bang time 
case. Even at low z, this is unacceptably high. 



E. More Complicated Profiles 

One may now ask whether changing the specific forms 
of £b(t*) and k(r) that we have used so far affects our 
results. We showed in Section [IV D| that very poor agree- 
ment with the kSZ data is obtained for the models that 
best-fit the supernova, CMB, and Ho data, but in Sec- 
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FIG. 6: The value of AT/T as a function of redshift for the 
best fitting models to the supernova, CMB and H data sets, 
from Section [IVB[ The corresponding value of f3 is shown on 
the right-hand axis. The median and 68% and 95% confi- 
dence intervals are shown as the black line, and the light gray 
and dark gray bands, respectively. The actual distribution is 
bimodal, and here we show only the models with low AT/T. 
Even for low redshifts, /3 is a large fraction of the speed of 
light. 
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FIG. 7: The 68% and 95% likelihood contours in the space of 
AT/T and Ho = lOO/i kms -1 Mpc -1 , for off-center observers 
at z = 0.05, 0.10, 0.15 and z = 0.20, for the MCMC sample 
constrained by SN+CMB+if data. 



The MCMC was constrained by the supernova, CMB, 
Ho, and kSZ data simultaneously. 

A plot of kSZ AT/T against Ho for these models is 
shown in Fig. [8j and may be compared with Fig. [7[ The 
models that maximize the likelihood have a (3(z) pro- 
file that is almost flat over the redshift range of interest 
(slightly larger at small z), and much less discrepant with 
the kSZ data. Relatively narrow spatial curvature and 
bang time profiles are preferred, extending out to only 
z ~ 0.2, and the bang time profiles are shifted towards 
the negative r direction (i.e. < in Eq. (|16|)). A 
preferred Hq of only 44.0 kms _1 Mpc _1 is obtained, and 
the fit to supernova and CMB data is also poor; the best- 
fit model is inconsistent with the data, with A% 2 ~ 270 
compared to ACDM. We conclude that this is because 
the fit is most sensitive to the kSZ data; it is easy to find 
models that are wildly inconsistent with the kSZ data, 
as evidenced by Fig. [7| and so models that minimize 
the x 2 with the kSZ data above all else are preferred. 
These tend to have low Hubble rates and narrow density 
profiles, features that are difficult to reconcile with the 
supernova and Ho data. As such, it seems that even with 
the significantly more complex bang time profile, a good 
fit to all of the data simultaneously is not possible. 

To further investigate the sensitivity of our results to 
the choice of profile parametrization, we now consider the 
model found by Celerier et al in [8] that was constructed 
to reproduce the ACDM values of luminosity distance 
and Hubble rate as a function of redshift, but without 
dark energy. The density profile on a hypersurface of 
constant t takes the form of a "hump" in this model, 
rather than a void, and the bang time gradient at low z 
is negative, so there are no shell crossings at early times. 
Instead, this model has double-valued redshifts at high 
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tion |IV C] we found that a good fit to the kSZ data could 
be obtained if a more complicated bang time profile was 
used. To see if a good fit to all of the observables is pos- 
sible with a more complicated model, we ran an MCMC 
simulation using the spatial curvature profile of Eq. d6| 
and the extended bang time profile given by Eq. (16). 



FIG. 8: The 68% and 95% likelihood contours in the space of 
AT/T and Ho — lOO/i kms -1 Mpc -1 , for off-center observers 
at z — 0.15, 0.20 and 0.55, for the MCMC sample constrained 
by SN+CMB+# +kSZ data. The AT/T are much lower 
than for the MCMC sample in Fig. |7| but the Hubble rate is 
too low to be considered consistent with observations. 
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FIG. 9: Upper panel: Bang time function, ts(r), of the LTB 
model found in 8 which reproduces the Hubble rate, H(z), 
and luminosity distance, cIl(z), of ACDM. Lower panel: Ve- 
locity with respect to frame in which the CMB is isotropic, 
f3(z), for the same LTB model. \/3(z)\ rapidly becomes very 
large (/3 ~ — c for z > 0.6), and produces a very poor fit to 
the kSZ data. 



z (see Section II B), which almost always result in large 
AT/T because Zi n and z out in Eq. 14 differ significantly. 
As such, this model is also strongly disfavored by current 
kSZ data (see Fig. [9}. 

In fact, for large fluctuations in the bang time function 
we expect that there will always be a large dipole seen 
by off-center observers at some range of redshifts, corre- 
sponding to lines of sight that pass near to regions with 
a non-zero bang time gradient [68 . It therefore appears 
that one cannot simultaneously fit the supernova, CMB, 
Hq and kSZ observations with a single LTB model unless 
one is prepared to violate one or more of the assumptions 
made in Section HI 



V. DISCUSSION 

In order to rigorously establish that A ^ 0, and 
that the concordance model of cosmology is correct, a 
model is required within which observations can be in- 
terpreted. The homogeneity and isotropy of the Uni- 
verse on large scales is often assumed, and most anal- 
yses are performed using the highly symmetric FLRW 
solutions of general relativity. These are not the only 
viable ways to model the Universe, however, and re- 
cent advances in observational cosmology allow us to 
empirically test alternatives rather than relying on as- 
sumed symmetries of space-time on the largest scales. In 
this paper we used the spherically- symmetric, dust-only 
Lemaitre-Tolman-Bondi (LTB) class of general relativis- 



ts cosmological models, in their full generality, to test 
the radial homogeneity of the Universe on large scales. 
In particular, we consider the magnitude of the kine- 
matic Sunyaev-Zel'dovich (kSZ) effect, which measures 
the dipole anisotropy of the CMB through a shift in 
the spectrum of CMB photons reflected from hot gas in 
clusters of galaxies. The kSZ effect is sensitive to large- 
scale inhomogeneity, as observers inhabiting an inhomo- 
geneous universe would generically expect to see large 
anisotropies on their CMB skies. 

In Section [TTJ we introduced the theoretical framework 
for specifying these models and calculating observables in 
them. We defined a simple parametrization of the LTB 
radial function k(r) that governs radial inhomogeneity 
at late times, and the function £#(r) that governs it at 
early times. This has been used to investigate under- 
dense 'voids' which have previously been shown to pro- 
duce good fits to the supernova data. We have discussed 
potential problems with models that are inhomogeneous 
at early times (i.e. that have non-constant including 
the potential for disruptions around our observed surface 
of last scattering. In models with a positive radial deriva- 
tive of the bang time, t' B > 0, it is found that shell cross- 
ing singularities form at early times, which pushes the 
surface of last scattering to later times. In models with 
negative bang time gradient, t' B < 0, regions with a neg- 
ative radial Hubble rate, H2 < 0, form and the distance- 
redshift relation, r(z), ceases to be monotonic. These 
features have observable effects that ultimately lead to 
predictions of a large CMB dipole anisotropy at low red- 
shifts. 

In Section IIIII we reviewed some of the observational 
constraints that can be imposed on void models with con- 
stant bang time. Three key sets of observables were con- 
sidered: The distance moduli of supernovae, the small- 
angle CMB power spectrum plus local Hubble rate, and 
upper limits on the magnitude of the kSZ effect for in- 
dividual clusters of galaxies. Voids can fit the supernova 
data easily, but are unable to fit recent measurements of 
the CMB and Ho simultaneously (they predict a value 
of Hq that is far too low). Similarly, voids which fit the 
supernova data predict a large CMB dipole at redshifts 
up to z ~ 1, which is inconsistent with current kSZ mea- 
surements. 

In Section IIVI we considered the effect on these con- 
straints of allowing the bang time to vary. This resulted 
in a significant increase in the freedom of the models, 
and allowed the supernovae and CMB+ffo data sets to 
be fit simultaneously, even with our simple parametriza- 
tion of the LTB radial functions. Models with small kSZ 
signals, consistent with the data, were also found, but 
these required more complex profiles and gave worse fits 
to the supernova data. We then proceeded to combine 
all of the observational constraints, and found that voids 
which are able to fit the supernovae, CMB and H predict 
an extremely large kSZ effect which is orders of magni- 
tude greater than the measured upper limits. A joint fit 
to the supernova, CMB, Hq, and kSZ data with a signif- 
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icantly more complicated bang time profile also failed to 
produce good agreement with the data. 

We also argued that any void model with a significant 
bang time inhomogeneity will produce a large kSZ effect 
at some redshift. Given that a varying bang time is nec- 
essary to resolve the low-#o problem, it seems that the 
combination of supernovae, CMB+iiZo and kSZ data is 
enough to effectively rule out LTB void models that at- 
tempt to explain cosmological data without dark energy, 
subject to the assumptions made in Section [I] This goes 
some way toward demonstrating the homogeneity of the 
Universe on large scales. 

We used the dipole approximation to calculate the 
magnitude of the kSZ effect in our models, but this only 
holds if the dipole term dominates the anisotropy of an 
off-center observer's sky [IIBES]. Otherwise, higher mul- 
tipoles become important, and the dipole approximation 
overestimates the kSZ effect. The dipole will dominate 
as long as most lines of sight on the observer's sky pass 
through the void, as will be the case if the observer is 
firmly inside it, for example. For the models in Section 
IV C| and IV D[ and the Celerier et at. model in Section 
IV E[ the inhomogeneity extends out to z > 2, and so 
the dipole approximation will always be a good one for 
observers at z < 0.6, where the kSZ data lie. The models 
with more complicated bang time profiles considered in 



Section IV E have much narrower inhomogeneities, how- 
ever, and so we would expect the dipole approximation to 
be worse. These models predict low kSZ magnitudes and 
are close to being consistent with the kSZ data, so any 
overestimate due to the dipole approximation will have 
little effect on their total x 2 7 which is anyway dominated 
by the poor fits to the supernova, CMB and Ho data. 
We therefore conclude that the dipole approximation is 
sufficient for our purposes. 

As we have tried to make clear throughout, our results 
are subject to several caveats that are summarized in 
Section [TJ The first, that we are exactly in the center of 
a perfectly spherically symmetric void, serves to simplify 
our calculations but is clearly unrealistic as it fails to 
take into account angular variations in, for example, the 
galaxy distribution. This could affect observables such 
as the dipole anisotropy of the CMB seen by off-center 
observers, potentially weakening the constraints we have 
derived using the kSZ effect. Considering ourselves as 
off-center observers [19] [52] and introducing linear per- 
turbations [30j [53] would produce more realistic models 
and allow more observational data to be used (e.g. the 
matter power spectrum) at the expense of a significant 
increase in complexity. A better understanding of linear 
perturbations would also go some way toward addressing 
our second caveat, that the formation of the last scatter- 
ing surface must be in an approximately-FLRW region. 
In general, one could expect features such as the cou- 
pling of scalar and tensor modes in LTB perturbations 
to produce secondary effects such as large B-mode polar- 
izations [30]. This type of effect is completely absent in 
linear perturbation theory about FLRW backgrounds. 



A particular limitation of LTB solutions as cosmolog- 
ical models is that they contain only dust, and cease to 
be applicable when radiation becomes important. If we 
want to approximate the Universe as an LTB model at 
late times, we must therefore match it to an appropriate 
solution containing radiation at early times. Solutions 
involving separate inhomogeneous matter and radiation 
fluids [31], spatially- varying physical quantities such as 
the photon-baryon ratio [14] , and scale-dependent initial 
power spectra [32 have been considered, and serve to give 
some idea of the extra freedom that might be obtained in 
more general models. Specifically, altering the location 
and properties of the surface of last scattering can have a 
profound effect on the predicted kSZ signal [14] and the 
observed CMB [31, 32 , and if one is prepared to consider 
this additional freedom, then our present results should 
not be expected to hold. 

Finally, let us consider LTB models in the context of 
general inhomogeneity. Rather than allowing space-time 
to be described by a single LTB metric, it has been sug- 
gested that the LTB geometry could be used as an effec- 
tive geometry to model the scale dependence of inhomo- 
geneity after some averaging procedure has been applied 
to the fine-grained structure of the actual inhomogeneous 
geometry of the real Universe [54] . This is a considerable 
departure from the situation we have been considering 
here. In particular, if other observers are able to con- 
struct similar effective spherically symmetric geometries 
about their own locations then we should no longer ex- 
pect distant clusters to see a large dipole in their CMB 
sky. This would completely relax the constraints that 
can be imposed from observations of the kSZ effect. 
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of the Universe at late times (and thus have any bear- 
ing on the low-i^o problem) would result in Compton 
^-distortions many times larger that can be reconciled 
with current observational constraints. 
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